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Abstract 



The nonlinear cr-model for disordered interacting electrons is studied in spa- 
tial dimensions d > 4. The critical behavior at the metal-insulator transition 
is determined exactly, and found to be that of a standard Landau-Ginzburg- 
Wilson <^ 4 -theory with the single-particle density of states as the order pa- 
rameter. All static exponents have their mean-field values, and the dynamical 
exponent z = 3. dn/dfi is critical with an exponent of 1/2, and the electrical 
conductivity vanishes with an exponent s = 1. The transition is qualitatively 
different from the one found in the same model in a 2 + e expansion. 

PACS numbers: 71.30. +h, 64.60.Ak 
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Despite much progress in our understanding of the metal-insulator transition (MIT) in 
recent years [0 , the search for a simple order parameter description of this phase transition 
in the spirit of a mean-field or Landau theory has so far proven futile. Early attempts at a 
mean-field theory of the Anderson transition of noninteracting electrons [[| failed because the 
most obvious simple order parameter (OP), viz. the single-particle density of states (DOS) 
at the Fermi level turned out to be uncritical 0. As a result, the Anderson transition can be 
described only in terms of a complicated functional OP H. At the Anderson-Mott transition 
of disordered interacting electrons, on the other hand, the DOS is generally believed to be 
critical, but it is not obvious how to construct an OP description in terms of it. Indeed, 
our understanding of this transition is largely based on a generalized matrix nonlinear a 
(NLcr)-model f|, for which no mean-field type fixed point (FP) is known, and which has 
been analyzed in terms of a small disorder expansion near d = 2. As a result, there is no 
simple description of the MIT analogous to, say, Weiss theory of ferromagnetism. This is all 
the more remarkable because of the technical similarity, first noted by Wegner [||], between 
the NLcr-model description of the MIT and that of Heisenberg ferromagnets. 

In this Letter we show that the NLcr-model for interacting disordered electrons || pos- 
sesses a saddle point solution which has all of the characteristic features of a Landau theory 
with the DOS as the OP. We further show that this solution corresponds to a renormal- 
ization group (RG) FP which is stable for d > 4. This establishes the exact (for d > 4) 
critical behavior at the MIT in this model, which can be summarized as follows: Let t be the 
dimensionless distance from the critical point at temperature T = 0, Q the energy distance 
from the Fermi level, and Q the DOS. Q vanishes according to 

Q(t,fi = 0)~t /3 , Q(t = 0,Q) ~tt p/vz . (la) 

For the critical exponents (3, is, and z, and for the exponents 7 and 77 characterizing the OP 
susceptibility, we find 

is = = 1/2, 7=1, 77 = 0, z = 3. (lb) 

All thermodynamic susceptibilities show the same critical behavior, 

X (t,n = 0)~t 1 '*, X (t = 0,fi)~fi 1/3 , (2) 

where x can stand for the density susceptibility dn/dfi, the specific heat coefficient 7 = 
lim T ^ Cy{t, T)/T, or the spin suceptibility Xs- As an argument of susceptibilities, fl denotes 
the external frequency, and fl and T can be used interchangeably in a scaling sense. The 
electrical conductivity a vanishes according to 

a(t, Q = 0) ~ t, a(t = 0,Q) ~ fl 2/3 , (3) 

so the conductivity exponent s — 1. 

In what follows we first derive Eqs.(l) by explicitly constructing the saddle point solution, 
and then using RG techniques to show that it is stable for d > 4. We then use scaling 
arguments to obtain additional information about susceptibilities, which leads to Eqs.(|2|),(0). 
We consider the matrix NLcr-model of Refs. i.e. a Gaussian field theory for a hermitian 
matrix field Q(x) with constraints (Q(x)) 2 = 1, with 11 the unit matrix, and trQ(x) = 0. 
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Q is a classical field comprising two fermionic fields. It carries two Matsubara frequency 
indices n, m and two replica indices a, (3 (quenched disorder has been incorporated by means 
of the replica trick). The matrix elements Q^(x) are in general spin quaternions, with the 
quarternion degrees of freedom describing the particle-hole and particle-particle channel, 
respectively. The action can be written 

S[Q,A] = 7^J dxtr (A(x)[g 2 (x) - 1]) +tr(d x g(x)) 2 + 2H Jdx tr(flQ( 

~\ E ff«[g(x)oQ(x)] u . (4) 

u=s,t,c 

Here G is a measure of the disorder, Q is a diagonal matrix whose elements are the external 
Matsubara frequencies u> n , and H is proportional to the free electron DOS. K SittC are cou- 
pling constants describing the electron-electron interaction in the particle-hole spin-singlet, 
particle- hole spin-triplet, and particle-particle or Cooper channel, respectively. K s < for 
repulsive interactions. [Q o Q] s ,t, c denotes a product in frequency space which is given ex- 
plicitly in Refs. Hp]]. Notice that we have enforced the constraint Q 2 = II by means of an 
auxiliary matrix field A(x). The constraint tr Q = will be enforced explicitly at every 
stage of the theory. 

The model, Eq.([|), is an effective one which was derived to capture the essence of the 
low-lying excitations of the system, i.e. the low-frequency large-wavelength behavior. The 
soft (i.e., diffusive) modes are given by correlation functions of the Q nm with nm < 0, while 
the DOS is determined by the average of Qnn It is therefore convenient to separate Q 

into blocks: 

Qnm = e(nm)Q&(x) + e(n)6(-m)g^(x) + 0(-n)0(m)(g+)^(x). (5) 

The conventional treatment of the NLu-model ||5],[l]] proceeds by integrating out A(x), using 
the constraint Q 2 = II to eliminate Q, and expanding the action in powers of q. Here we use 
a different approach inspired by treatments of the 0(N)-NL<r-model in the large-N limit ||. 

Since the further development will be closely analogous to the treatment of an O(N)- 
Heisenberg model, let us pause to point out the similarities between these models. The 
matrix elements Q nm correspond to the massive a-component of the O(N)- vector field, while 
the q nm correspond to the massless 7r-fields. The disorder G plays the role of the temperature 
in the magnetic model, it is the control parameter for the phase transition. HQ is in some 
respects analogous to the magnetic field conjugate to the order parameter a. The last term in 
Eq. (f|) has no analogy in the Heisenberg model. We will see that in the mean- field treatment 
of the matrix model presented here it plays a rather trivial, although crucial, role. 

We proceed by integrating out the massless g-field. Since the action, Eq.(||), is quadratic 
in Q and hence in q, this can be done exactly. We then look for a saddle point of the 
resulting effective action S[Q, A]. This task is simplified by restricting ourselves to saddle 
point solutions that are spatially constant, diagonal matrices with diagonal elements Q", 
A". This ansatz is motivated by the fact that < Qnm > has these properties. For simplicity 
we restrict ourselves to the particle-hole spin-singlet channel, i.e. we put K t = K c = 0. We 
will see later that this restriction does not influence the critical behavior. We further use a 
short-range model interaction, so that K s is simply a number. Again it can be shown that 
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a Coulomb interaction leads to the same critical behavior 0. We find for the saddle point 
solution, 



(«) 2 = i + tE E 2k1VK " 



4 Vmtlib 2 + |(A^ + A«)] 5 



1 + n j^ 1 2~TGK., 



^ r»2 _|_ ifAa i A« ^| 
rti=0 r ' 2 \ "i ' ll n 1 -n+m) 



(6a) 



A* = 2GHu n /Q%. (6b) 

We discuss several aspects of this result. Firstly, Q"=o = Q, which is the DOS normalized 
by the free electron result, decreases with increasing disorder (remember K s < 0). This is 
the well known 'Coulomb anomaly' in the DOS and it is proportional to K s as well as 
G. Technically, this is due to the replica structure of the theory: All terms on the r.h.s. of 
Eq.(^ap that are independent of K s vanish in the replica limit and have not been shown. 
Physically, it reflects that fact that the Coulomb anomaly is due to the interplay between 
interactions and disorder. Secondly, Eqs.(j|) constitute an integral equation for Q% which 
involves an integration over all frequencies. As noted above, the NLcr-model is designed to 
describe only low-frequency behavior, and can not be trusted at high frequencies. However, 
on physical grounds it is clear that Q" tends to a constant at large frequencies, so for fixed 
K s Q will vanish at a critical value G c of G. We can easily determine the critical behavior 
for 5G = G - G c < 0: We find Q ~ (-5G) 1 / 2 , which is the first relation in Eq.([Ta|). 
Furthermore, the leading low- frequency behavior for G = G c involves only integrals up to 
the external frequency, which is the region where the theory is controlled. For d > 4 we find 
Q(Cl) ~ fi 1 / 3 , which is the second relation in Eq. (|Ta|) with the exponent values given in Eq. 

We have shown that our saddle point actually corresponds to a minimum of the free 
energy by expanding the action to second order in the fluctuations 8Q,SA about the saddle 
point. Details of this calculation will be reported elsewhere 0. The result is a positive 
definite Gaussian matrix, so the saddle point is indeed a minimum. One can then integrate 
out 5 A to obtain the order parameter correlation function G nin2 n2n4 (k). At zero frequency, 
and close to the critical point (i.e., for Q — >• 0) the result simplifies substantially and we find 



G/4 
k 2 + Q^/a 



Goo,oo(k) — u2 ■ ^ 2 /c2 , (7) 



where the bare correlation length £ is given in terms of a complicated integral which is finite 
for d > 4. From Eq.(|7]) we obtain three more critical exponents: u = (3 = 1/2, rj = 0, and 
7 = 2u = 1, cf. Eq.([lt]). 

Apart from the DOS we are interested in the transport properties. Let us first consider 
the charge diffusion coefficient D c . As a hydrodynamic quantity it can be obtained from 
the NLcx-model by a direct calculation of the particle-hole spin-singlet c/-propagator [HQ. 
Inserting the saddle point solution, Eqs.@, into Eq.(||), one reads off the q- vertex function, 
and a matrix inversion yields the corresponding propagator. The latter has the usual dif- 
fusion pole structure, and D c is obtained as the coefficient of the momentum squared. We 
find 

Dc= G(H + K S Q) • (8) 
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We see that at the transition D c vanishes like the OP. The behavior of the conductivity can 
be obtained by multiplying D c with dn/dfx. The latter is less straightforward to obtain, 
since as a thermodynamic quantity it is not simply given by the hydrodynamic behavior. 
We will determine its critical behavior from Eqs. flTTD below. 

We now use RG techniques to show that these results represent the exact critical behavior 
of the model in d > 4. The RG will also provide an easier route to a derivation of Eqs.@ and 
(0) than a direct calculation would be. Let us return to the action, Eq.(^j). If we integrate 
out q, we get Eq.(f|) with Q replaced by Q plus terms obtained by contracting q- fields. Of 
the latter, one term is linear in A. The resulting action can be written in the form 

S[Q,A] = ~cj dx(<9 x Q(x)) 2 - fdx tr(tA(x)) + 2H J dx *r(fiQ(x)) 

7rT r 
—T £ tf„[Q(x).Q(x)] u + Ul / dx tr(A(x)Q 2 (x)) 

4 u=s,t,c 

+ / (ixj^ A/(x)(m2)/jAj(x) + (other terms). (9) 
i,J 

Here c — u\ — 1/2G and t is a matrix composed of — (1/2(2) H. and the term of 0(A) coming 
from the g-contractions. I = (nm, a/3, i) where i labels the spin-quaternion space, and u-i is 
a matrix which is finite for d > 4. The explicit expressions for t and 112 will not be needed 
for our present purposes, they will be given elsewhere 0. The 'other terms' in Eq. (|j|) all 
come from contracting g-fields. They can easily be constructed diagrammatically, but will 
turn out to be irrelevant for our purpose. 

We now apply standard power counting to the action, Eq.(0). Our parameter space 
is spanned by \x = {c, t, H, K s>tjC , u\, U2}, and the coupling constants of the 'other terms.' 
In looking for a RG fixed point (FP), we follow Refs. || and || in fixing the exact scale 
dimensions of our fields Q and A to be [Q] = d/2 — 1, [A] = d — 2 (we define the scale 
dimension of a length to be -1). This corresponds to fixing the exponents r] and fj defined 
by the wavenumber dependence of the Q- and A-correlation functions to be rj = and 
fj — d — 4, respectively. With [Q] = [T] = d, we find the scale dimensions of the various 
coupling in Eq.® to be [c] = 0, [t] = 2, [H] = 1 - d/2, [K sM ] = 2 - d, [m] = [u 2 ] = 4 - d. A 
power counting analysis of all of the 'other terms' in Eq.(||) is straightforward. The result 
is that the coupling constants of all of these terms have scale dimensions which are smaller 
than 4 — d for d > 4. We conclude that the Gaussian FP given by fi* = {c, 0,0,0, ...} is 
stable for d > 4. The only relevant parameter is t, and the correlation length exponent v is 
v — l/[t] = 1/2. c is marginal, as expected, and all other parameters, including H and K s ^ c ,i 
are irrelevant (the combinations HQ and K s t c T are, of course, relevant, which reflects the 
fact that a finite frequency or wavenumber drives the system away from the critical point). 
This Gaussian FP obviously corresponds to the saddle point solution discussed above. It 
also is, structurally, just the ordinary Gaussian FP of 4 -theory: For scaling purposes, A is 
equivalent to Q 2 . This is a consequence of our choice of the scale dimensions of Q and A, 
and becomes explicit if one neglects the irrelevant 'other terms' in Eq. (^[) and integrates out 
A. The stability of the FP proves that the critical behavior obtained from the saddle point 
solution is exact in d > 4. 

The RG arguements given above imply that the order parameter obeys the scaling rela- 
tion 
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Q(t, HQ, u h u 2 , ...) = h x - d l 2 Q{iti tv ', HQb 1+d / 2 , Ml 6 4 - d , u 2 b 4 ~ d } ...), (10) 

where v = 1/2, and we have suppressed all parameters with scale dimensions smaller than 4- 
d. The exponents j3 and z follow from Eq. (|10D by using standard arguments [f| . The crucial 
point is that u\ and u 2 are dangerous irrelevant variables: Solving the theory explicitly in 

1/2 

the saddle point approximation, we see that Q(t, ft — 0, U\ — > 0) ~ u x , and Q(t = 
0, — ► 0,u 2 — > 0) ~ {u 2 /u^) 1 ^. Equation ([10]) then immediately yields (3 = 1/2, z = 3 in 
agreement with Eq. ([Tb]) . We emphasize that K s is noi dangerously irrelevant, even though 
Q(t = 0,Q) is proportional to \[K~ S . The point is that the vanishing of the (bare) K s just 
shifts the transition point to infinity, making the mean-field transition inaccessible, but does 
not change the critical behavior. 

We finally derive Eqs. (0), @. To this end we notice that all susceptibilities have been 
identified in terms of the coupling constants of the NLa- model: dn/dfi ~ H + K S , \ s ~ H + 
K t , 7 ~ H |T0|,|I|j. H, K s , and K t are all irrelevant, which means that all three susceptibilities 
vanish at the critical point. Furthermore, K s t scale to zero faster than H, so the asymptotic 
scaling behavior of all three susceptibilities is the same and given by that of H, or 7. The 
only remaining difficulty is to correctly incorporate the dangerous irrelevant variables. Let 
us consider the singular part f s of the free energy density /, which satisfies the scaling 
equation 0, 

f a (t,T, Ul ,...) =b-^f s (tb 1 /»,Tb z ,u 1 b 4 - d ,...). (11a) 
In the critical regime we have f s ~ Q 2 ~ l/«i, and hence 

f s (t,T)=b-^f s (tb^,Tb z ). (lib) 



We see that hyperscaling breaks down in the usual way: d in Eq. ( 11a ) gets replaced by 4 



By differentiating Eq. ( [[ 1 b| ) twice with respect to temperature, and using z = 3, we find 

7 (t,T) = r 1 7 (t6 2 ,T6 3 ), (12) 

and hence Eqs.(Q) [|1^,|T3|]. Finally, combination of this result with Eq.(|j) yields Eq.(||). 
We note that the criticality of dn/dfi is very remarkable, since it does not show in the 
2 + e expansion treatment of the same model |],[IJ] , while it is an important feature of the 



Mott-Hubbard transition, which recently has enjoyed a revival of interest |TT . 

Several questions arise from these results: (1) Can the technical approach used here 
be extended to d < 4? An obvious possibility is to introduce an (artificial) parameter iV 
multiplying the soft modes, and perform a 1 /iV-expansion in analogy to the 0(iV)-model 
0. (2) What is the connection between the present results and those obtained in d = 2 + e, 
where dn/dfi is found to be uncritical in perturbation theory ||14|| ? A possible answer is 



that dn/dfx has an essential singularity near d = 2 which is invisible in the 2 + e expansion. 
Another possibility is that the Gaussian FP discussed here is not continuously connected to 
the one found near d = 2. (3) Can the present techniques be extended to the more general 
field theory [0] underlying the NLa-model? If that is the case, one would expect the upper 
critical dimension of the resulting OP description to be 6, or, for technical reasons similar 
to those in Ref. @, 8. Which of the two models would capture more of the physics relevant 
in d = 3 would a priori be unclear. These problems will be pursued in the future @ . 
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In conclusion, we have shown that the generalized NLa-model for interacting disordered 
electrons, which has been studied extensively in d = 2 + e, contains a metal-insulator tran- 
sition with mean-field like critical behavior for d > 4. This transition shows qualitatively 
different features from the one obtained in d = 2 + e: It is order parameter driven rather 
than soft-mode driven, and it shows a critical dn/d/i. 

This work was supported by the NSF under grant numbers DMR-92- 17496 and DMR- 
92-09879. 
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